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ABSTRACT 

We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and 
angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass 
loss rate drops by a factor of wl.5-2, while the angular momentum loss rate drops by a factor of 
«4 as the distance decreases past the Alfven surface. This reduction in angular momentum loss is 
due to the interaction between the stellar and planetary Alfven surfaces, which modifies the global 
structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the 
opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum 
loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The 
interaction however, begins at greater separation if the overall field topology of the star and the planet 
are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars 
harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss 
can compete with, and perhaps dominate, spin-up due to tidal interaction. 
Subject headings: stars: magnetic field — stars: coronae — planet-star interactions 



1. INTRODUCTION 

Nearly two decades after the first discovery of a planet 
orbiting a main-sequence star outside of the solar system 
(Mayor & Queloz 1995), hundreds of other exoplanets 
have been observed (Mayor et al. 2003). Many of these 
planets have masses in the Jupiter class, and orbit their 
host star at distances of <10 stellar radii. Such massive 
close-in planets are also called "hot Jupiters" . 

Planets that are far out in the stellar system cannot 
significantly affect their host star via tidal (gravitational) 
interaction. They also cannot affect the star via magnetic 
interaction (even if the planet has a strong internal mag- 
netic field) since the super- Alfvenic stellar wind prevents 
the propagation of information back from the planet to 
the star along the field lines; the star is thus essentially 
blind to the existence of the planetary magnetic field. 
The conceptual boundary which separates the sub- and 
super- Alfvenic regimes is called the "Alfven surface" lo- 
cated at a radial distance from the star known as the 
" Alfven radius" . This surface is defined as the collection 
of critical points at which the wind flow speed, u, is equal 
to the Alfven speed, ua, or in other words, the Alfvenic 
Mach number, Ma — u/ua = 1. In the case of close-in 
planetary orbits at distances comparable to or less than 
the stellar Alfven radius, magnetic interactions are more 
easily facilitated and can, at least in principle, induce a 
feedback on the star via star-planet magnetic interaction. 

Signatures of star-planet interaction have been ob- 
served through local modulations in the Ca II K emission 
line — an indicator of chromospheric activity — in several 
planetary systems (Shkolnik et al. 2005a, b, 2008, see also 
Hartman (2010)). In addition, an apparent enhancement 
in stellar activity in stars with close-in planets compared 
to stars with planets located at greater orbital radii ap- 
pears in the form of an overall increase in observed X-ray 
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(Kashyap et al. 2008, see, however, Poppenhaeger et al. 
2010) and UV (Shkolnik 2010) flux. 

Recently, Pont (2009) performed an empirical study 
between stellar rotation rate, orbital semimajor axis, 
stellar and planetary radii, and the planet-to-star mass 
ratio. This study suggests that stars that harbor hot 
Jupiters rotate faster than stars with planets located at 
larger orbital separation, or stars without planets at all. 
Pont (2009) followed (Mardling & Lin 2002) in suggest- 
ing that the excess in stellar angular momentum is due 
to tidal torques that synchronize the stellar rotation with 
the planetary orbital motion. It is also possible to view 
this result as a reduction in the angular momentum lost 
from stars with hot Jupiters, rather than as the host stars 
being spun up. (Lanza 2010) suggests that in many ob- 
served planetary systems the synchronization time-scale 
is too long for tidal interaction to cause this excess in 
stellar angular momentum. Instead, he proposes that 
the angular momentum excess is due to star-planet mag- 
netic interaction that leads to a reduction of magnetic 
breaking and angular momentum loss to the stellar wind 
(Weber & Davis 1967). He concludes that both mecha- 
nisms can take part in the process of angular momentum 
transfer and that the dominant of the two mechanisms 
depends on the particular set of stellar-planetary param- 
eters, as well as the stellar type. 

In this report, we present MHD numerical simulations 
of the star-planet interaction, where we focus on calcu- 
lating the stellar angular momentum loss rate for a rep- 
resentative set of stellar and planetary parameters as a 
function of the planetary orbital radius. We describe the 
simulation details in §2, describe the results and in §3. 
We summarize in §4. 

2. SIMULATION 

The simulations described here were performed in the 
same manner as in (Cohen et al. 2009b). The model used 
solves the set of MHD equations, based on the surface dis- 
tribution of the stellar magnetic field, and provides a self- 
consistent, steady-state solution for the stellar corona 
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and stellar wind (Cohen et al. 2007). The planet is rep- 
resented by a dipole magnetic field with boundary con- 
ditions typical of Jupiter. Since the simulation provides 
the solution for the density, velocity, magnetic field, and 
pressure, we can calculate the realistic Alfven surface and 
hence the stellar angular momentum loss rate (see Cohen 
et al. 2009a for a detailed description of this procedure). 

We adopt some of the stellar properties of the well ob- 
served HD 189733 system (#* = O.76i? , M* = O.82M 
(Mayor et al. 2003), and P rot = 11.95c? (Fares et al. 
2010)). (Fares et al. 2010) have recently investigated the 
surface magnetic field distribution of HD 189733; here we 
adopt an idealized dipolar stellar magnetic field in order 
not to complicate any modulation of model results with 
planetary distance. The planetary field is assumed to be 
dipolar as well, with equatorial field of B p = 1 G. The 
choice of a field somewhat lower than that of Jupiter is 
based on the expected slower rotation of close-in plan- 
ets as a result of tidal locking (Sanchez-Lavega 2004), 
and the observed inverse correlation between planetary 
magnetic field strength and rotation period (e.g. Durand- 
Manterola 2009). 

In addition to the stellar Alfven surface, we expect 
the planet to possess an analogous surface due to plane- 
tary outflow. In this simulation, we mimic the planetary 
outflow by introducing the planet as a mass source for 
the plasma outflow, where the flow rate depends on the 
particular choice of the planetary boundary conditions. 
The planetary radius used here is R p = 1.95 Rj (Rj is 
Jupiter's radius) with boundary conditions for the den- 
sity, n v = 10 7 cm~ 3 , and temperature, T p — 10 4 K. 
This particular choice of parameters produces a plane- 
tary outflow of about 50 km s^ 1 . 

Here, we study modulations in stellar angular momen- 
tum loss rate as a function of the semimajor axis, a. 
Once the stellar Alfven surface is defined, we calculate 
the mass and angular momentum loss rates which depend 
on the topology of this surface. We explore values of a 
ranging from 12R* to 3i?*, and calculate the loss rates 
for each case. We omit the case of a = 5R±, since at 
this distance the planet is located exactly at the stellar 
Alfven surface. This causes the stellar Alfven surface to 
be an open surface with the grid space of the simulations 
and the integration required for the loss rate calculations 
is ill-conditioned. We also test how the orientation of the 
planetary magnetic field affects the loss rates by defin- 
ing a planetary dipole magnetic field which is aligned or 
anti-aligned to the stellar dipolar magnetic field. 

3. RESULTS & DISCUSSION 

The MHD solution for the particular case of an aligned 
planetary dipole and semimajor axis a = 12R* is illus- 
trated in the left panel of Figure 1. Also shown with an 
arrow is the range of semimajor axis investigated. The 
six panels on the right of Figure 1 show the stellar and 
planetary Alfven surfaces for different values of a and for 
aligned and anti-aligned dipoles. 

The results for the mass and angular momentum loss 
rates as a function of a are shown in Figure 2. It can be 
seen that the mass loss rate drops by about a factor of 
1.5-2 at a = 9R± for the aligned case and at a = 11-R* 
for the anti-aligned case. The trends in the plot for the 
angular momentum loss rate are similar, but here the 
change reaches a factor of 4 at a — 8R±. The reason for 



this change is clearly seen in Figure 3, where we show 
the mass flux distribution, as well as the location of the 
Alfven surface over the star-planet meridional plain for 
the cases of a = 12R* and a = 6R±. In the case of 
a = 12R+, the front of the planetary Alfven surface is 
beyond the stellar one, so that the existence of the planet 
does not affect the structure of the stellar corona. In the 
case of a = 6R±, the planet is located inside the stel- 
lar Alfven surface, so the stellar corona and stellar wind 
the "feel" the obstacle (i.e. the planet and its magneto- 
sphere) , and the large-scale structure of the corona is dis- 
rupted. In particular, both on the planetary side and on 
the opposite side, the mass flux near the equatorial region 
is significantly reduced due to the existence of the planet. 
This explains the 50% reduction in mass flux when the 
planet affects the corona. The fact that the coronal side 
opposite to the planet is modified demonstrates that the 
effect of the planet on the corona is global. We would 
like to point out that another reference simulation with 
different set of stellar and planetary parameters resulted 
in the same trend, so that this effect does not depend on 
the particular choice of parameter set. In addition, an 
initial result from a simulation of the HD 189733 plan- 
etary system, which includes realistic stellar magnetic 
field, as well as a planetary orbital motion resulted in a 
total mass flux of 6T0~ 14 M Q j/r _1 , which is close to the 
value of 4 • 10~ 14 M© yr" 1 obtained from this simulation 
for a = 9R* (as observed for HD189733). 

Since in our case the star and the planet are tidally 
locked, the planet blocks a significant amount of the 
stellar wind, and prevents the stellar wind from open- 
ing up field lines and escaping (Cohen et al. 2009b). For 
the particular set of parameters used in our simulations, 
the planetary magnetosphere has a large enough cross- 
section to block the stellar wind over a significant area. A 
planet without a significant internal magnetic field would 
not only be too small to affect the mass flux through the 
Alfven surface, but would probably suffer from a strong 
erosion by the stellar wind. In our simulations, the radius 
of the cross section of the planetary magnetosphere is not 
larger than about 1-2 stellar radii. Therefore, the area 
which is blocked by the planet is not more than about 
25-30% of the total area of the stellar Alfven surface. 
Nevertheless, the whole coronal hemisphere that faces 
the planet is affected, so that the stellar mass loss rate 
decreases by approximately a factor of 1.5-2 as a result 
of the interaction with the planetary magnetosphere. 

The four right panels of Figure 2 show that the planet 
begins to affect the corona when the planetary Alfven 
surface, (and not the planet itself) starts to interact with 
the stellar one. They also shows that the two surfaces 
merge at greater separation for the case of an anti- 
aligned planetary dipole as compared to the aligned case. 
This explains the smaller value of a at which a drop in 
mass and angular momentum loss rates compared to the 
aligned dipole case is seen. This merging of the surfaces 
at smaller a is due to particular magnetic field lines that 
connect the planetary magnetic field and the stellar one. 

The amount of change in the loss rates depends mostly 
on the size of the planetary magnetosphere, which is de- 
termined by the particular choice of planetary param- 
eters. The alignment of the planetary field has only a 
minor effect on the change in angular momentum loss. 
This is due to the fact that the wind is blocked primarily 
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by the magnetic pressure of the planetary field, which 
depends on the magnitude of the field but not on its po- 
larity. It is most likely that a tilted planetary dipole will 
produce only slightly different result, depending on the 
tilt angle. 

Pont (2009) presented empirical evidence for "excess" 
rotation in stars hosting close-in planets. This was inter- 
preted in terms of tidal forces acting towards spin-orbit 
synchronization and transfer of orbital to rotational an- 
gular momentum. While such angular momentum ex- 
change doubtless takes place, the MHD models presented 
here demonstrate that the stellar angular momentum loss 
rate through magnetized winds also likely plays a signif- 
icant role. Close-in planets reduce this angular momen- 
tum loss and reduce stellar spin-down. The effect is large 
enough — a factor of up to 4 — that it could dominate the 
observed rotation excesses. As noted by Pont (2009), 
spin-up by tidal effects requires substantial orbital decay 
since the dissipation of the protoplanetary disk. The re- 
duced magnetic braking effect found here would imply 
that much less severe orbital decay would be required. 

Full understanding of the distribution of planetary or- 
bits and stellar rotation rates likely requires considera- 
tion of both tidal and magnetic angular momentum loss 
and exchange. Improvements in the MHD modeling pre- 
sented can be made at the expense of substantial compu- 
tational cost by including time-dependent effects of the 
star-planet interaction using more realistic surface mag- 
netic field distributions. 
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4. CONCLUSIONS 

In this report we calculate variations in stellar angu- 
lar momentum loss rate of a star hosting a hot Jupiter 
planet as a function of the planetary orbital distance. We 
show that once the stellar and planetary Alfven surfaces 
interact with each other, the stellar wind topology in the 
hemisphere facing the planet changes and the angular 
momentum loss to the wind decreases. The simulation 
also shows that this interaction begins at a greater dis- 
tance if the overall field topology of the star and the 
planet are opposite, though the magnitude of change in 
angular momentum loss rate depends mostly on the size 
of the planetary magnetosphere. This reduction in angu- 
lar momentum loss for close-in planets likely contributes 
to, and perhaps even dominates, the observed excess ro- 
tation found in stars with close-in planets. 
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Fig. 1. — Left: Steady-state solution for the stellar corona and the planet for the a = 12R*, aligned dipole case. The star and the planet 
are shown as yellow and blue spherical shades, respectively. Red lines represent the stellar magnetic field lines, while yellow lines represent 
the planetary magnetic field lines. Background color contours are of number density and the Alfven surfaces of the two bodies are shown 
as solid white lines. The pink arrow marks the direction at which the planetary location has been modified with distance. Right: Stellar 
and planetary Alfven surfaces for the cases of a = 12, 11, 8i?* (left to right, respectively) colored with contours of the local mass flux at 
each point. Top row shows the solution for the aligned planetary dipole case while the bottom row shows the solution for the anti-aligned 
planetary dipole case. 
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Fig. 2. — Mass loss rates, Mdot> (top, units of solar mass per year) and angular momentum loss rates, Jdot, (bottom, units of [g cm 2 s 2 ]) 
for the case of aligned planetary dipole (solid blue line, triangles) and anti-aligncd planetary dipole (dashed red line, squares). 
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Fig. 3. — Mass flux distribution (color contours) and the location of the Alfven surface (black line) displayed on the y=0 (star- planet) 
meridional plain. The solution for the a = 12R* case is shown on the top panels while the solution for the a = 6-R* case is shown on the 
bottom panels. Left panels show the global solution with the star as a red sphere, and the planet as a blue sphere (a = 12i?*) or yellow 
sphere (a = 6-R*). The right panels show a zoom on the vicinity of the stellar corona hemisphere which is not occupied by the planet. In 
both right panels, the Alfven surface for the a = 12R* case is shown as a solid line, while the Alfven surface for the a = 6R* case is shown 
as a dashed line. 



